Astronomy & Astrophysics manuscript no. 
(will be inserted by hand later) 



Spin Period Evolution of Recycled Pulsar in Accreting Binary 

J. Wang 1 , CM. Zhang 1 , Y,H, Zhao 1 , Y. Kojima 2 , H. X. Yin 3 , and L.M. Song 4 



o 

(N 

> 
O 

CO 
CN 

x 

Oh- 

6 

CO 



1 National Astronomical Observatories, Chinese Academy of Sciences, Beijing 100012, China, 
e-mail: zhangcm@bao .ac.cn 

2 Department of Physics, Hiroshima University, Higashi-Hiroshima 739-8526, Japan 

3 School of Space Science and Physics, Shandong University, Weihai 264209, China 

4 Institute of High Energy Physics, Chinese Academy of Sciences, Beijing 100049, P. R. China 
the date of receipt and acceptance should be inserted later 

Abstract. We investigate the spin-period evolutions of recycled pulsars in binary accreting systems. Taking both the accretion 
induced field decay and spin-up into consideration, we calculate their spin-period evolutions influenced by the initial magnetic- 
field strengths, initial spin-periods and accretion rates, respectively. The results indicate that the minimum spin-period (or 
maximum spin frequency) of millisecond pulsar (MSP) is independent of the initial conditions and accretion rate when the 
neutron star (NS) accretes ~> 0.2M o . The accretion torque with the fastness parameter and gravitational wave (GW) radiation 
torque may be responsible for the formation of the minimum spin-period (maximum spin frequency). The fastest spin frequency 
(716 Hz) of MSP can be inferred to associate with a critical fastness parameter about a> c = 0.55. Furthermore, the comparisons 
with the observational data are presented in the field-period (B - P) diagram. 
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1. Introduction 

Neutron stars (NSs) are usually detected as either normal pul- 
sars (single or in a binary) with magnetic-field B ~ 10 12 G and 
spin-period P ~ 0.5 s, or millisecond pulsars (MSPs) (half in 
binary) with B ~ 10 8 5 G and P ~ 20 ms (e.g. Bhattacharya 
& van den Heuvel 1991; Lorimer 2008). The spin-periods and 
surface magnetic-fields assigned to these two different systems 
are found to span different ranges of B - P values (Hobbs & 
Manchester 2004; Manchester et al. 2005; Lorimer 2008), and 
the updated statistical distributions of B and P are plotted in 
Fig.Q](data from ATNF pulsar catalogue). In this figure, we no- 
tice that the distributions of both B and P are almost bimodal, 
as firstly suggested by Camilo et al. (1994), with a dichotomy 
between "normal" pulsars (P ~ 0.1 s - 10s and B ~ 10 11 G 
- 10 13 G) and MSPs (P ~ 1 ms - 20 ms and B ~ 10 8 G - 
10 9 G). This abnormal phenomenon has led to the idea that ra- 
dio pulsars in binary systems have been "recycled", i.e. they 
have been spun up due to mass accretion during the phase of 
mass exchange in binaries (Alpar et al. 1982; Radhakrishnan & 
Srinivasan 1982; Bhattacharya & Srinivasan 1995). The weak 
magnetic-field of MSP strongly supports the idea that it decays 
in the binary accretion phase (e.g. Taam & van den Heuvel 
1986; van den Heuvel 2004). The currently widely accepted 
view about this phenomenon is the accretion-induced evolu- 
tion, i.e., NSs accreting materials from their low-mass binary 



companions are spun up by the angular momentum carried by 
the accreted materials, while the magnetic -field decays (Taam 
& van den Heuvel 1986; Bhattacharya & van den Heuvel 1991; 
van den Heuvel 2004). The X-ray binaries are the evolutionary 
precursors to these "recycled" MSPs. It is evident that B and 
P of X-Ray pulsars and recycled pulsars are correlated with 
the duration of both the accretion phase and the total amount of 
matter accreted (Taam & van den Heuvel 1986; Shibazaki et al. 
1989). If the NS accretes a small quantity of mass from its com- 
panion, e.g. ~ 0.001M o - 0.01M o , a recycled pulsar with the 
mildly weak field and short spin-period (B ~ 10 10 G, P ~ 50 
ms) will be formed (e.g. Francischelli, Wijers & Brown 2002), 
like PSR 1913+16 and PSR J0737-3039 (Lyne et al. 2004). 
The direct evidence for this recycled idea has been found in 
low mass X-ray binary (LMXB) with the accretion millisec- 
ond X-ray pulsar (AMXP), e.g. Sax J 1808.4-3658 (Wijnands 
& Klis 1998), and in observing the transition link from an X- 
ray binary to a radio pulsar PSR J1023+0038 (Archibald et al. 
2009). 

At the end of the accretion phase of NS/LMXB (accreted 
mass ~> O.2M ), the NS magnetic-field may arrive at a bottom 
value about 10 8 G and spin-period may reach a minimum about 
millisecond (Cheng & Zhang 1998, 2000), remaining a MSP 
(van den Heuvel & Bitzaraki 1995ab). Since then, spin-down 
due to dipole radiation has been conceived as the "standard evo- 
lution" of MSP (Urpin & Konenkov 1997; Urpin, Geppert & 
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Fig. 1. Distribution of magnetic-field and spin-period for 1864 pulsars (data from the ATNF pulsar catalogue, see Manchester et 
al. 2005). Here, 141 binary pulsars and 13 magnetars (Kouveliotou 2003) are separately identified for the statistics of magnetic- 
field and spin-period. The 24 accretion millisecond X-ray pulsars discovered by RXTE (see e.g. Wijnands 2005; Lamb & Yu 
2005; Yin et al. 2007) are put in the spin-period histogram. 



Konenkov 1997, 1998). The switch-on of a radio pulsar could 
exist, whose radiation pressure is capable of ejecting the ac- 
creted matter out of the system and prevent any further accre- 
tion. This 'radio ejection' phase in the binary will account for 
the formation of a MSP (Burderi et al. 2002ab). However, the 
question whether all observed MSPs could be produced within 
this recycled framework has not been quantitatively addressed 
until now. The accretion induced collapse (AIC) of a white 
dwarf as an alternative mechanism for the formation of MSP 
has been studied recently (Ferrario & Wickramasinghe 2007 ab; 
Hurley et al. 2009; Wickramasinghe et al. 2009). Recently, an- 
other idea for MSP formation is proposed, which claims that 
these objects are born with low magnetic fields (Halpern & 
Gotthelf2010). 

Then the mechanisms for the accretion-driven field decay 
and spin-up in binary NSs have been suggested by many re- 
searchers in various proposals. The accretion flow and thermal 
effects speed-up Ohmic dissipation of NS crust currents, which 
accounts for the field decay and spin evolution (Geppert & 
Urpin 1994; Geppert, Page & Zannias 1999; Urpin & Geppert 
1995; Romani 1990). The interactions between the accreted 
matter and magnetic-field in the course of accretion phase are 
studied to account for the field decay (Melatos & Phinney 
2001; Payne & Melatos 2004; Konar & Bhattcharya 1999; 
Konar & Choudhury 2004; Lovelace, Romanova & Bisnovatyi- 
Kogan 2005). Not only is the total amount of accreted mass 
considered to influence the final field strength and minimum 
spin-period (van den Heuvel 1995; van den Heuvel & Bitzaraki 
1994, 1995ab), but also the influences of accretion rate may be 
a factor (Wijers 1997; Cumming et al. 2001; Cumming 2005; 
Zhang & Kojima 2006). 

Taking the accretion-induced field decay into considera- 
tion, we investigate the formation of the high spin frequency 
(short spin-period) for a MSP, and the influences on it by the 
initial conditions (e.g. initial spin-period and initial field) and 
accretion rate. The spin-up torque and field evolution make 
their spin-period and field distributions in B - P diagram to 



change from the regions of B ~ 10 13 G and P ~ 1 s - 100 
s to that of B ~ 10 8 ~ 9 G and P ~ 1 ms - 20 ms as observed. 
This paper is organized as follows. In Section 2, we review the 
main equations that dominate the B-P relation in the accretion 
phase. The effects by the initial conditions and accretion rate on 
the spin-period are presented in Section 3. Section 4 contains 
the summary and discussion. 

2. Field Versus Period Relation 

2.1. Description for the Model 

As declaimed by the accretion induced magnetic field decay 
model (Zhang & Kojima 2006), if the magnetic field is suf- 
ficiently strong, e.g. ~ 10 12 G, and the star spin very slowly 
initially, the accreted matter will be channeled onto the polar 
patches by the field lines, where the compressed accreted mat- 
ter causes the expansion of magnetic polar zone in two direc- 
tions, downward and equatorward, which makes the magnetic 
flux in the polar zone diluted. With accretion, the field decays 
and magnetosphere shrinks until the magnetosphere reaches 
the surface of NS, where a bottom field is obtained to be about 
~ 10 s G. 

Based on the above model, the accretion-induced field evo- 
lution is obtained analytically with the initial field B(t — 0) = 
B (Zhang & Kojima 2006), 



(1 - [C expi-y) - 1] 2 )I 

Here, we define the parameters as follows, y — j^f-, the 
accreted mass AM = Mt, the crust mass M cr ~ O.2M , 
C = 1 + ^1 - xl ~ 2 with x\ = (|f) 4/7 . B f is the bottom 
magnetic-field which is defined by the NS magnetosphere ra- 
dius matching the stellar radius, i.e., Rm(B/) = R. Rm is defined 
as Rm - <PRa where Alfven radius Ra (Eisner & Lamb 1977; 
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Ghosh & Lamb 1977) reads, 

R A = 3.2 x 10 8 (cm)M- 2/7 //3o 7 n^ 1/7 , (2) 

The model dependent parameter <p is about 0.5 (Ghosh & Lamb 
1979b hereafter GL; Shapiro & Teukolsky 1983; Frank et al. 
2002). Myj is the accretion rate in units of 10 17 g/s. /^o is the 
magnetic moment in units of 10 3() Gem 3 . The mass m = M/M Q 
is in the unit of solar mass. Using the relation RM(Bf) = R, we 
can obtain the bottom field, 

Bf = 1.32 x 10 8 (G)(-^-)Mfl 6 (3) 
Mis 

where M l% = M/lO l& g/s and R 6 = fl/10 6 cm. As declaimed 
by Inogamov & Sunyaev (1999), in the accretion process, a 
boundary layer forms between the innermost disk and magne- 
tosphere due to the transition for rotating velocity of plasma 
from Keplerian to the spin velocity of NS. When the accretion 
rate increases, this layer may spread (Inogamov & Sunyaev 
1999), and as a result more matter is accreted to the polar 
gap and diffuses to the entire surface of the NS. This effect 
will influence on the field decay efficiency of the model by 
Zhang & Kojima (2004), since some spreading matter has 
little contribution to the field lines dragging of the polar 
patches. Moreover, such a spreading will expand the polar 
cap area, which makes the polar cap to occupy the entire 
NS surface while the magnetosphere does not reach the NS 
radius. The bottom field of NS is determined by the condi- 
tion that the magnetosphere equals the NS radius (Zhang 
& Kojima 2004). If there is no field decay while the magne- 
tosphere radius is bigger than NS radius, then the modified 
bottom field by considering the spreading of the polar cap 
will be bigger than the one obtained for the ideal frozen 
plasma in polar cap zone. Namely, the magnetic field of re- 
cycled NS is slight bigger than the ideal value of the model 
by Zhang & Kojima (2004). 

To study the spin evolution of NS in accretion phase, we 
employ the formula for accretion-induced spin-up given by GL, 

-P = 5.8X10- 5 [A-**?T£] 

M 6 

xBl^PLlM^) syT 1 , (4) 

where we define the parameters, the surface field B\i = 
B/lO l2 G, the moment of inertia I45 = 7/10 gem 2 , the X-ray 
brightness (L = GMM/R) L37 in units of 10 31 erg /s, respec- 
tively. The dimensionless parameter n(u> s ) is the fastness pa- 
rameter which is explained in the following subsection. 

2.2. Fastness Parameter 

For a slowly rotating magnetic NS, the rotation of the star cou- 
ples to the orbital motion of the disk matter at the Alfven radius 
when accretion flow is approximately radial. So the limit of an- 
gular velocity (Q s = 2n/P) for a NS can reach the Keplerian 

velocity (Q^CRm) = ). In general, the spin velocity 

is always less than the orbital velocity, i.e. Q, <k GI^Rm) = 
Therefore, the relative importance of stellar rotation 



can be described by the ratio parameter of the angular veloci- 
ties (Eisner & Lamb 1977; Ghosh & Lamb 1977; Li & Wang 
1996, 1999; Shapiro & Teukolsky 1983), 

oj s = _^L_ = i.35[(— r 2/ X 5/7 ]fi^p-%- 3/7 , (5) 

which plays a significant role in our entire understanding of 
accretion to the rotating magnetic NSs. 

To calculate the accretion torque which acts on a magnetic 
NS accreting matter from a disk, GL introduced a dimension- 
less accretion torque n(u> s ) depending primarily on the fastness 
parameter. A simple expression for n(u> s ) is given by GL, 

n(co s ) = 1.4 x(i^^), (6) 

where (o c is the critical value. For a slowly rotating star (co s <K 
1), GL found that the dimensionless function n(co s ) decreases 
with lo s and become negative for o> s > io c if n(oi s ) ~ 1 .4. In 
addition, GL also stressed (±> c ~ 0.35 for their model. However, 
the subsequent work (Ghosh & Lamb 1992) indicates that co c 
is unlikely to be less than 0.2, but it could be as large as 0.9. 
By the aid of this dimensionless parameter, we can establish 
the relation between the total torque on the star and the torque 
communicated to the star by the magnetic-field lines that thread 
the inner transition zone. 

In the following subsections, we study the spin evolution 
of NS in a binary system, where the influence by fastness is 
included. 

3. Spin Evolution of Accreting NS 

In this section, we investigate the spin-period evolution of ac- 
creting NSs by solving equations ([T), ©, © and (O with dif- 
ferent initial conditions and accretion rates, where we set the 
usual parameters for NS, e.g. m = 1.4, R^ = 1.5, 4> - 0.5 and 
Mcr = 0.2M G . 

3.1. Influence by Initial Conditions 

We consider the spin evolution for a wide range of initial condi- 
tions (initial magnetic-fields and initial spin-periods): (1). vary- 
ing the initial magnetic-field 5 x 10 11 G, 5 x 10 12 G and 5 X 10 13 
G while setting spin-period and accretion rate as 100 s and My], 
respectively (Fig. [2j; (2). varying the initial spin-period 1 s, 10 
s & 100 s with a certain field (5 x 10 12 G) and constant accretion 
rate (M n ) (Fig.0). 

For the samples, we classify them as normal pulsars, MSPs 
and binaries. From the figures, we can see that most MSPs are 
consistent with the binary pulsars, while only a few binaries fall 
in the same region of the normal pulsars that have a magnetic- 
field ~> 10 11 G. Recycled pulsars are expected to be found 
only to the right of the spin-up line and to the left of the "death- 
line" (Bhattacharya & van den Heuvel 1991), which is shown 
indeed the case in these figures (Fig. 2, Fig. 3) for almost all 
the binaries and MSPs. This strongly supports the view that the 
MSPs have been recycled (Alpar et al. 1982; Radhakrishnan 
& Srinivasan 1982). Firstly, with a small amount of mass trans- 
ferred, the NS is spun-up from the death valley with long period 
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Fig. 2. The diagram of the magnetic-field versus spin-period 
for pulsars. The spin-up line represents the minimum spin- 
period to which such a spin-up may proceed in an Eddington- 
limited accretion, while the "death-line" corresponds to a po- 
lar cap voltage below which the pulsar activity is likely to 
switch off (Bhattacharya & van den Heuvel 1991). The evo- 
lutionary tracks in B - P diagram are plotted with different ini- 
tial magnetic-field strengths, but with same initial spin-period 
(P = 100 s) and same accretion rate M\i\ The solid, dash and 
dot lines denote the numerical solutions for equations ([TJ, ©, 
© and © with initial field B„ = 5x 10 13 G, 5 xlO 12 G and 5 
xlO 11 G, respectively. 
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evolves to the hundred millisecond regime, and here the system 
accretes ~ 0.0 1M . After the system accretes about ~ 0.2 M , 
the field reaches ~ 10 8 G - a bottom value, and the spin period 
enters the millisecond regime, which is inferred corresponding 
to the condition ai s — w c . Namely, the bottom spin period or 
maximum spin frequency can be obtained by 

VsMax = V k (R M = R)0Jc , (7) 

where v k (R M = R) = 1300(Hz)(M/lAM e ) l/2 (R/15km)- V2 is 
the maximum spin frequency at stellar surface. If the critical 
fastness parameter is u> c = 0.55, then the maximum spin fre- 
quency can be as high as 715 Hz, which is similar to the fastest 
spin frequency of MSP observed (716 Hz). The lower the criti- 
cal fastness, the lower the final spin frequency achieved. 

For our numerical calculation, we take the effect of fastness 
parameter into consideration. If the critical fastness parameter 
is set as 0.9, the numerical solutions for equations ([TJ, ©, (0 
and d6) suggest that the evolutionary curves in the B—P diagram 
are a little below the spin-period equilibrium line. However, 
the effect of the fastness parameter is significant when the B-P 
evolutionary track is close to the equilibrium period line (spin- 
up line). This indicates that the spin angular velocity of the star 
tends to reach the Keplerian angular velocity at the inner edge 
of the accretion disk, but never exceeds this angular velocity. In 
fact, the primary effect of the fastness parameter is to force the 
evolution curves back to or below the equilibrium period line 
and ensure that the evolutionary track cannot go beyond this 
line. 

From these tracks, we can see that despite of the differences 
for initial fields and initial spin periods, the NSs can approach 
at a certain minimum spin period and bottom field after accret- 
ing ~> 0.2M o . If the accreting phase is over, a recycled pulsar 
can be found in the region of B - P diagram with B ~ 10 8 " 9 G 
and P ~ 10 ms, forming a MSP. On account of the existence of 
minimum spin-period and bottom field (e. g. B ~ 10 8 G), the 
P-distribution and B -distribution imply that the MSPs accumu- 
late close to the region of low field and short period, which is 
an explanation of the bimodal distributions for all pulsars, as 
shown in FigQ] In short, for LMXBs the different initials give 
the similar final value of spin-period, which is insensitive to the 
initial conditions. 



Fig. 3. The same meaning as Fig. [2] with different initial spin- 
periods (P - 1 s - dot line, 10 s- dash line, 100 s - solid line) 
at the fixed initial magnetic-field (B„ — 5 x 10 12 G) and same 
accretion rate Mn. 

and strong field, e.g. the NS in HMXBs (B ~ 10 12 G, P ~ 10 
s) like Her X-l (P = 1.24 s, B = 3 X 10 12 G, see van der Meer 
et al. 2007). With the long-lived accretion phase the binary ac- 
cepts sufficient mass from its companion and the NS yield a 
substantial field decay, as in the case of MSPs (Bhattacharya 
& van den Heuvel 1991), e.g., the fastest known pulsar PSR 
J1748-2446 (B ~ 10 8 G, P = 1.4 ms, see Hessels et al. 2006). 

As shown in Fig. 2, three B-P tracks with the different ini- 
tial field values begin to follow one track when the spin period 



3.2. Influence by Accretion Rate 

In this part, we study the influence by accretion rates on the 
spin-period evolution. We set the accretion rate as M\^,M\i 
and Mi g at the fixed magnetic-field (5 x 10 12 G) and spin-period 
(100 s), respectively(see Fig. |4j». From the Fig.|U we find that, 
with different accretion rate, the spin evolution follows differ- 
ent tracks. If we set the critical fastness parameter as 0.9, the 
tracks with Eddington accretion rate (M\% = 1) matches the 
evolution of most recycled pulsars well, but the track with M^, 
lies below the "death-line". Moreover, the curves with a high 
accretion rate and high critical fastness have the tendency to 
approach the spin-up line. 

For a comparison of the magnetic-field and spin-period 
distributions between theory and observation, we plot the 
statistical distributions of magnetic-field and spin-period 
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Fig. 4. The same meaning as Fig. [2] with different accretion 
rates, at the fixed initial magnetic-field (B — 5 x 10 12 G) and 
initial spin-period (P„ = 100s), from the upper to bottom, M^ 
- solid line, M\q - Dash line and Mif, - dot line, correspond- 
ing to the bottom fields of ~ 10 7 5 G, ~ 10 8 G and ~ 10 85 G, 
respectively. The critical fastness parameter is set to be (±> c - 
0.9. 



obtained from the calculation for theoretical model with 
various initial conditions (see caption of Fig. |5). 100 sam- 
ples of final B-field and final spin period (B, P ) have been 
done, which are plotted to compare the observed data of bi- 
nary pulsars as shown in Fig.|5J where we notice a bimodal 
distributions of magnetic-field and spin-period of recycled 
pulsars (initial distributions of B and P at right parts of 
plots compared to the final evolved values at the left parts 
of plots). As compared to the observational distributions of 
binary pulsars, we notice that the both histograms for the 
MSPs (left parts of Fig. |5), theoretical and observed data, 
are very similar. For the spin period, more theoretical data 
than the observed values occurs at the regimes close to mil- 
lisecond, which may imply that more MSPs may not evolve 
to their minimum spin period of about millisecond. 

3.3. Spin-period Versus Accretion Mass (P - AM) 

To investigate the spin-period evolution with accretion mass, 
we consider the influences by initial conditions and accretion 
rate and plot the evolutionary curves in P - AM diagrams (see 
Fig. 01. As a detailed illustration, Fig. [6] implies the follow- 
ing results. (1). After accreting certain mass (e.g. ~ 0.01M o ), 
the spin-period evolution seems independent of the initial 
magnetic-field. (2). The influence on P— AM by the initial spin- 
period on the spin evolution exists when AM <~ 0.001M G and 
disappears when the accreted mass AM > 0.00 1M . When a 
small amount of mass (~< 0.001M Q ) is transferred, the spin- 
period mildly changes (Fig.|6]l, which may yield a HMXB with 
the spin-period of some seconds, such as Her X-l (van der 
Meer et al. 2007) and Vela X-l (Quaintrell et al. 2003). After 
a conveniently long time, with the weight of material accreted 



~ O.O1M , NS is spun-up to a shorter period (P ~< 100ms) 
which corresponds to the system of DNS, e.g. PSR 1913+16 
(Breton 2009). The binaries with longer accretion phase, e.g. 
LMXBs, will accrete sufficient mass (AM ~> O.1M0) from 
their companions and yield a lower magnetic -field and shorter 
spin-period (P ~< 20 ms), as in the case of MSPs, e.g. SAX J 
1808.4-3658 (Wijnands & van der Klis 1998; Wijnands 2005) 
and PSR J 1748-2446 (Lorimer, 2008). (3). At the beginning 
of accretion, NSs with different accretion rates follow differ- 
ent curves in the P - AM plane. After accreting ~ 0.01M o , the 
spin-period approximately follow the same track with different 
accretion rate, until to the same minimum period. Therefore, 
NSs in LMXBs, for the Z source with Eddington luminosity 
and Atoll source with less luminosity, e.g. ~ 1% Eddington ac- 
cretion rate (Hasinger & van der Klis 1989), should have the 
similar spin frequencies while accreting ~ O.2M . 



4. Summary and Discussion 

Considering the influences exerted by the different initial 
magnetic-field, initial spin-period and the accretion rate, we in- 
vestigate the evolutionary tracks of recycled pulsar in B - P 
diagram. The main results are listed below. 

(1) . From the evolutionary tracks in B - P diagram, we find 
that the magnetic-field of NS decays with accreting material, 
and the spin-up evolution proceeds at the same time. The pri- 
mary effect of the fastness parameter is to force the evolution- 
ary curves to reside below the equilibrium period line. The crit- 
ical fastness (oj c ) can infer the final spin-period of recycled NS 
(maximum spin frequency), for instance, u c = 0.5 will corre- 
spond to a spin of 650 Hz (1.54 ms), which is near the highest 
observed value of MSP and AMXP 

(2) . The influences of the initial magnetic-field, spin-period 
and accretion rate are tested, respectively. The evolutionary 
curves present that the spin-period (frequency) decreases (in- 
creases) with the accretion and the minimum spin-period is in- 
sensitive to these two initial conditions when it is spun-up to 
P ~ 10 ms (v ~ 100Hz) after NS accreting about ~> O.2M . 

(3) . It is found from the B-P diagram, a few observed MSPs 
locate above the "spin-up line", while the trajectories based 
on our model always lie below this line. Firstly, owing to the 
Ohmic dissipation, the buried field of MSP may re-emerge after 
accretion, which can account for a slightly increase of field of 
the star (e.g. Young & Chanmugam 1995; Bhattacharya 2008). 
Secondly, the magnetic model we use here is based on the 
"dilution of polar magnetic flux" due to accretion, which is 
idealized. Due to numerous plasma instabilities, the Rayleigh- 
Taylor instability at the disk-magnetosphere interface because 
of the high-density disc matter being supported against by the 
low-density magnetosphere plasma and Kelvin-Helmholtz in- 
stability because of the discontinuity in the angular velocity of 
matter at the boundary (Ghosh & Lamb 1979a), these instabili- 
ties may result in the penetration of the magnetosphere, prying 
the field lines aside and azimuthally wrapping the field lines by 
the disc matter (e.g. Romanova 2008; Kulkarni & Romanova et 
al. 2008), which may modify the field strength evolution then 
perturb the spin evolution. 
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Fig. 5. Plots of histograms for spin-period and magnetic-field obtained from the calculation for the theoretical model and for the 
observed binary pulsars (red-shaded). In theoretical calculations, we assume that the initial B-field follows a Gaussian distribu- 
tion, centered at 5 x 10 12 G, and ranged from 5 x 10 10 G (e.g. Hartman et al. 1997) to 10 14 G (e.g. Kaspi 2010); As for the initial 
spin period, it is assumed to be Gaussian distribution, centered at 30 s, ranged from 0.5 s to 100 s; The different accretion rates 
are set from to Mig. The theoretically evolved B-field and spin period values occurred at the left parts of plots, which are 
similar to the observed distributions. 
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Fig. 6. Plots of spin-period versus accretion mass with different conditions, i.e., initial magnetic-fields (left), initial spin-periods 
(middle) and accretion rates (right), where the critical fastness parameter is set to be a> c = 0.9 in calculation. It is noticed that the 
initial conditions (Bo and Pq) and accretion rate has little influence on the P - AM plot after the system accretes ~> 0.01M o ; In 
other words, the spin periods of NS/LMXBs should be primarily related to how much mass accreted. 



(4). The formation of the fastest spin frequency of MSP is 
an interesting topic (Lorimer 2008; Chakrabarty et al. 2003; 
Chakrabarty 2005). The known fastest radio pulsar is PSR 
J1748-2446ad (Ter 5) at 716 Hz (Hessels et al. 2006), whereas 
the highest AMXP spin frequency in LMXBs is 619 Hz (e.g. 
Wijnands 2005). However, the recent report of 1122 Hz burst 
oscillation frequency, interpreted as a spin frequency, in XTE 
J 1739-285 (Kaaret et al. 2007) has been reported but not yet 
confirmed on account of the statistical significance (actual sig- 
nificance is only 3<x). Now the detected spin is less than the 
believed breaking spin frequency of pulsar of about ~ 1000 Hz 
(Lattimer & Prakash 2004), then what causes the special max- 
imum spin frequency of MSP? 

The spin evolution is affected by the accretion torque, 
which consists of two parts (Ghosh 1995; Frank et al. 2002): 
the stresses associated with matter accreting from the inner 



edge of the disk (A^o = M(GMRm)^) and that associated with 
the magnetic field coupling the star with the disk. The total 
torque N can be conveniently expressed in terms of A^o and 
the dimensionless torque. The spin-down torques which bal- 
ances the spin-up torque will also contribute to forming the 
maximum spin frequency. Slow rotators (u> s < cl> c ) are spun 
up (« > 0) by this accretion torque. o) s < 1 is presumed, 
since for co s > 1 the steady accretion is of course impossible. 
Here, (±> c of approximately ~ 0.3 - 0.5 is the critical fastness 
at which the torque changes direction, thus the maximum spin 
frequency v s - 2jtQ. s - a) c Vk ~ (0.3 — 0.5)v£, where vu = 2nD.it. 
The observed maximum upper kHz QPO frequency is 1330 
Hz (van der Klis 2000; Zhang et al. 2006) that is interpreted 
as the Keplerian frequency vk of inner disk orbit, so this im- 
plies a maximum spin frequency to be about ~ 400 - 665 Hz 
which is consistent with observations. To limit the occurrence 
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of high spin of MSP over 1000 Hz, the braking torques have 
been considered. Electromagnetic braking torque can brake the 
stellar spin (Ghosh 1995, 2006), which are usually expressed 
in the form n oc ju 2 (Q,/c) 3 , where n is the magnetic dipole 
moment of the star. In NS/LMXB, the mass quadrupole can 
occur on account of the accretion that arises the gravitational 
wave (GW) radiation (Bildsten 1998; Melatos 2007; Melatos 
& Payne 2005; Melatos & Peralta 2010; Ghosh 2006; Vigelius, 
Payne & Melatos 2008; Watts et al. 2009). The angular mo- 
mentum is lost from the emission of GW which may spin-down 
the NS. With accretion, the gravitational wave and accretion 
torques balance, but the total torque is effectively negative be- 
cause the gravitational wave torque increases with accreting. 
Hence, there might arise a maximum spin frequency of several 
hundred Hertz for a MSP. 
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